On Geometric Properties of Passive Random Advection 
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We study geometric properties of a random Gaussian short-time correlated velocity field by 
considering statistics of a passively advected metric tensor. That describes universal properties of 
fluctuations of tensor objects frozen into the fluid and passively advected by it. The problem of 
one-point statistics of co- and contravariant tensors is solved exactly, provided the advected fields do 
not reach dissipative scales, which would break the symmetry of the problem. Asymptotic (t — > oo) 
duality of the problem is established, which in the three-dimensional case relates the probabilities 
of the volume deformations into "tubes" and into "sheets" . 



I. INTRODUCTION 

A problem of passive advection in a turbulent medium 
attracts considerable attention as a solvable model of tur- 
bulence. Exact solutions can be found in a simplified 
case, when the velocity field is chosen to be a random, 
short-time correlated Gaussian process. Statistics of den- 
sity, concentration, passive vectors advected by such a 
field were investigated by many authors (see, e. g., 
where intermittent nature of the fluctuations, non-trivial 
scalings of structure functions and anomalous role of the 
dissipation were discovered. All these features are very 
common in the general picture of turbulence and, there- 
fore, the problem of passive advection can serve as a 
model for developing corresponding analytical tools. 

In the present paper we consider passive advection (in 
the Lie sense) of a second-rank covariant tensor in d- 
dimensional space. Though our master equation for the 
probability density function (PDF) of the tensor (Eq. (^) 
below) is very general, we concentrate mainly on statis- 
tics of a symmetric (metric) tensor gij . One-point statis- 
tics of any tensor object frozen into the fluid can be re- 
lated to statistics of such a tensor. We do not impose 
any restrictions (such as incompressibility) on the veloc- 
ity field, and, therefore, statistics in both Eulerian and 
Lagrangian frames are studied. Also, we are only inter- 
ested in the "initial stage" of the advection, when the 
advected field does not reach dissipative scales. This al- 
lows us to explore the symmetries of the problem, which 
are broken when dissipation is included. 

We show that the probability-density function of the 
eigenvalues of the metric is governed by a d-particle 
Hamiltonian that can be split into two non-interacting 
parts. Its non-universal part describes the motion of the 
center of mass (the determinant g of the metric) and can 
be separated from the motion relative to the center of 
mass, i. e. dynamics of the metric's eigenvalues normal- 
ized to their geometrical mean, Xi/g 1/d . The Hamilto- 



nian of the latter motion is of the Calogero-Sutherland 
type, remains the same in both Lagrangian and Eulerian 
frames of reference, and therefore describes the universal 
properties of the advection. These properties are dictated 
by the symmetry of the problem. The exact integrability 
of the Calogero-Sutherland Hamiltonian is known to be 
related to SL(d) symmetry: the Hamiltonian can be rep- 
resented as a quadratic polynomial in terms of the gen- 
erators of the corresponding algebra j^^^l- The eigen- 
functions of this Hamiltonian are the so-called Jack poly- 
nomials, which are symmetric homogeneous functions of 
the eigenvalues. This allows us to find exactly all mo- 
ments (T m ) of any tensor T advected by the fluid. In- 
deed, calculating any such moment reduces to averag- 
ing expressions of the type Tr k (g n ), which are symmetric 
polynomials in terms of the metric's eigenvalues, and can 
therefore be expanded in Jack polynomials of degree nk. 
We illustrate this method by calculating exactly all mo- 
ments of passively advected vectors and covectors, in par- 
ticular, of the magnetic field in kinematic regime and of 
the passive-scalar gradient. We also demonstrate how 
this approach works in the general case of a passively 
advected tensor of any rank. 

Calculating the moments requires knowing the statis- 
tics of the metric g with special initial conditions, gij{t = 
0) = 5ij. However, it is also interesting to consider 
the evolution of the PDF of the symmetric tensor g^ 
subject to arbitrary initial conditions. In this context, 
we show that a beautiful dual picture exists: the time- 
dependent PDF of the tensor becomes asymptotically 
(t — > oo) invariant under the inversion of the eigenvalues 
with respect to their geometrical mean. For example, in 
three dimensions, that means that if a magnetic field ad- 
vected by ideally conducting fluid develops flux tubes, it 
must develop magnetic sheets with the same probability. 

The paper is organized as follows. In Section ||, we 
derive the master equation for the PDF of the met- 
ric's eigenvalues, and analyze the symmetry properties of 
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this PDF. In Section II], we present a simple method of 



transforming the PDF between Eulerian and Lagrangian 
frames, which is important in the case of a compress- 
ible velocity field. Section [j^ discusses general properties 
of solutions for the PDF in two- and three-dimensional 
cases. In Section |v|, we show how the symmetry of the 
problem allows to calculate all the moments of passively 
advected tensors. The paper is written in a self-contained 
manner, all the necessary definitions and derivations are 
summarized in the Appendices. 



II. MASTER EQUATION 

A covariant second-rank tensor field ipij (t, x) passively 
advected by the velocity field £ fe (i,x) evolves according 
to the following equation: 



d t <pij + £ k <Pij,k + t k t <Pk 3 + Z*<pik = 0, 



(1) 



where = d£ k /dx % , and <Pij,k — d<pij/dx k . Let £*(t,x) 
be a Markovian Gaussian field: 

(e(*,x)^(t',x'))=^(x-x')5(t-t'), 
/s« (y) ~ KoS" - k 2 (yW + 2ayV ) , y - 0, (2) 

where a is the compressibility parameter, and «2 = 1 f° r 
simplicity. Here a can vary between — l/(d+ 1) for the 
incompressible flow and 1 for the fully compressible flow. 

In order to determine the statistics of the tensor, we 
follow a standard procedure P,|l0| and introduce the 
characteristic function of ip(t, xj: 



Z(t,a) = (exp{cr y v3 4i (t,x)}) . 



(3) 



This function is a Fourier transform of the PDF of the 
matrix elements tfij . Clearly, Z is independent of x due 
to spacial homogeneity. We find that Z satisfies 



d t Z 



[1 + a(d + 1)] a ij 



r<'J . 



d 



2 V ,)cr> 



dZ 
da> k 



2aa lJ - 



d 
d 



, dZ 
da mn 

d 



da k 



da 



Ik 



d_ 

da 1 



Jk 



da u 



Z, 



(4) 



where d is the dimensionality of space. This equation was 
derived by taking the time derivative of Z, using Eq. (|l]), 
and splitting Gaussian averages. We obtain the equa- 
tion for the probability density function of (p by Fourier- 
transforming (0): 



P(<p) = J exp {-tff*V«} z i?) n dfjr ' 



(5) 



The original equation (|l]) preserves symmetry proper- 
ties of the tensor ipij , which means that we may restrict 



our consideration either to advection of symmetric or an- 
tisymmetric tensors. Both reductions can be done in a 
similar fashion. For our present purposes we only con- 
sider fluctuations of a symmetric covariant tensor. The 
corresponding results for a contravariant tensor are sum- 
marized in Appendix |c[ We will use both (co- and con- 
travariant) pictures when discussing statistics of passive 
vectors in Section [v| 

In the symmetric case, the PDF (|) can be factorized 
as follows: 



(6) 



where g is the symmetric part of (p. One may think of 
the tensor g as of a metric associated with the medium. 
Due to spacial isotropy, P depends only on the eigenval- 
ues Ai , . . . , Xd of g. After rather cumbersome but essen- 
tially simple calculations, we establish the following mas- 
ter equation for the PDF of the eigenvalues of the metric: 

5 t P = 2(2a + l)^A?^+2a^AiA ' W> 



dX 2 1 — Z^'^QXidXj 
dP 



[3d + 4 + 2a(d 2 + 3d + 3)] ^ A, — — 



dX t 

, (dP dP 

+{a+1) ^—Ad*rdx- 

+ ^d(d+l)(d + 2)[l + a(d+l)]P (7) 

(from here on the overtildes are dropped). Among the so- 
lutions of this equation, those corresponding to the PDF 
must be non-negative, finite, and normalizable. The nor- 
malization is as follows llll : 



j dAi . . . dA d P(Ax, . . . , \ d ) JJ \Xi - A,-| = 1 



(8) 



Clearly, the original stochastic equation ([!]) preserves the 
signature of the metric. We will restrict ourselves to the 
case of all positive A's. Since there is no means of distin- 
guishing between different orderings of the eigenvalues, 
the PDF must be a symmetric function with respect to 
all permutations of Ai, . . . , X ( i- 

We should now notice that in logarithmic variables 
Zi = log(Aj), the master equation (Q) describes the dy- 
namics of d pair-wise interacting particles on the line. 
Furthermore, we can consider these dynamics in the 
reference frame associated with the center of mass of 
the particles z = z i- Denoting the coordinates of 
the particles in this frame Q = Zi — z, and noticing 
that det(g) = g = exp(zd), we find that P now satis- 
fies 



d t P = d[l + a(d + l)] <^ 2g 



, 8 2 P 
dg 2 



(2d + 5) g 



dP 
dg 



2 



+ -{d + l)(d + 2)P} +2(1 + a) 



d 



d 2 P 



^ dc 2 + 2 ^ 

i=l i<] 



1 



dP dP 



(9) 



where the <i variables Ci)--->Cd are n °t independent. 

d =0. The Hamiltonian remaining after the dynam- 
ics of the center of mass are separated, is translation- 
ally invariant, therefore the total momentum of the par- 
ticles 'YlidP/dQ) is conserved. The normalization rule 



dCi-dC d 6(Ci + --- + ( d )\J(0\ / dgg 



Jd-l)/2p_ 



AO- 



2<i(d-l)/2 



d 1 

| J sinh - (Q 



p = l, 



(10) 



where by C we denote the set {0, . . . , Cd}- The operator 
in the square brackets is a Sutherland Hamiltonian H$, 
which is exactly solvable (see, e.g., ^^2 L|; this Hamil- 
tonian appeared in a similar context in []5 j) . The Hamil- 
tonian Hs is the same for co- and contravariant tensors, 
and in both Eulerian and Lagrangian frames. 

It is important that Hs is self-adjoint with respect to 
the measure (fLOf). Its eigenfunctions are the so-called 
Jack polynomials, that are homogeneous polynomials 
in exp(£j) and are symmetric with respect to all per- 
mutations of Q. Their construction is discussed in Ap- 
pendix We will use particular eigenfunctions of this 
operator in Sec. [v|. 

We see that if P is initially chosen in a factorized form, 
P = Pi(g)P2(Ci, ■ • ■ , Cd), it will remain so factorized at 
all times. Thus, the statistics of g are independent of 
the statistics of the £'s at all times if they are initially 
independent. In particular, this property of Eq. (^) al- 
lows us to consider separately the PDFs for the deter- 
minant of the metric and for the logarithmic quanti- 
tiesC l = log(A l5 - 1 /'i): 



(11) 



An additional symmetry emerges in this context: (0) 
and ( |Io|) remain invariant if the coordinates Zi of all par- 
ticles are simultaneously reflected with respect to their 
center of mass. Such reflection leaves the center of mass 
intact and reverses the signs of all i. e. transforms 
all \i into g 2 l d j\. The origin of this symmetry can be 
understood if we notice that the master equations for 



S(t,g) = J d d (\J(()\P(gX), 
F(t,Q= [ dgg^/ 2 P(g,()- 



the PDFs of the dimensionless quantities Git 



-l/d 



and G 



ik 



gi/ d gik are samej although the initial 



9ik 



stochastic equations are different. This symmetry leads 
to nontrivial results for d > 3, and will be considered 



in Section IV 



III. EULERIAN AND LAGRANGIAN PDF'S 



The equation for the metric-determinant PDF S(t, g) 
follows from Eq. (||): 



,d 2 S 
dg 2 



dS 1 



d t S = 2g 2 —j + (2d + 5)g^ + ^-(d + l)(d + 2)S, (12) 



where we have rescaled time by the factor of 7 = d [1 + 
a(d+l)]. This factor is always non-negative and vanishes 
if the velocity field is incompressible, a = — l/(d+ 1), in 
which case any time-independent function S(g) is a solu- 
tion. Note that the right-hand side of Eq. (fL2h becomes a 
full derivative when multiplied by the Jacobian g( d ~ 1 ')/ 2 . 
The solution of this equation is a log- normal distribution: 



S(t,g) 
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-(rf+l)/2 



y/E^jt 



■ exp 



(\og(g)+ 1 ty 
8jt 



(13) 



where we took the initial distribution in the form 
S(0,g) = %-l). 

This result can be simply understood if we note that 
the determinant g obeys the same equation as p 2 , the 
squared density of the medium. The density satisfies the 
continuity equation, which can be written in logarith- 
mic form: 



5 t logp + ^5 fc logp + e fe fe = 0. 



(14) 



Since the time increments of £ fc are independent iden- 
tically distributed random variables, the Central Limit 
Theorem implies the normal distribution of log p. In- 
deed, either from Eq. ( |l2| ) or directly from Eq. (fil|), 
one can easily establish that the density PDF R(t, p) — 
2p d S{t,p 2 ) satisfies d t R = (7/2) (p 2 R)". 

So far, we have worked in the Eulerian frame, con- 
sidering statistics at an arbitrary fixed point x. Now 
we show how the one-point joint Eulerian and La- 
grangian PDF's are related. Let us assume that ini- 
tially Lagrangian particles are uniformly distributed in 
space. We denote the Eulerian PDF Pe(p, Ci ^ x )> tne 
Lagrangian PDF Pl(p, (', t, y), where y is the Lagrangian 
label (initial coordinate of the Lagrangian particle), 
and p = I det(<9y/<9x)| (the density of the medium). The 
relation between Pe and Pl can be established from the 
following: 

P E (p, C; t, x) = (S(p- p(t, x)) S (( - C(t, x))> = 

^(6 (x - x(t, y)) 6 ( P p(t, y)) S (( - C(t, y))>. (15) 
P 

Since the one-point PDF Pe(p, Ci^ x ) is independent 
of position (due to spacial homogeneity), we can inte- 
grate (|l5|) with respect to x. Also noting that the one- 
point PDF Pl(p, Ci y) i s independent of y, we get: 
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p e (pX) = -Pl(pX)- 
p 



(16) 



Transformation to the Lagrangian frame can also be per- 
formed on the level of the original stochastic equations 
such as (H) with the aid of the stochastic calculus (see, 

e.g., gg). 

In our considerations, if we choose intially 5(0, g) cx 
5(g — pg), we may substitute p = ^fg in formula (jig). We 
see therefore that only the PDF of g is affected by the 
transformation between Eulerian and Lagrangian frames. 
The Lagrangian version of S(g) is: 



S(t,g) 



-(d+l)/2 



■ exp 



(iog(g)-7*r 

8jt 



(17) 



Analogous results for the contravariant case are presented 
in Appendix |^. 

The log-normal statistics such as (13) and ( jl7j ) are a 
signature of this problem, and they will also be present 
for fluctuations of the eigenvalue ratios in asymptotically- 
free regimes, i. e. where different ratios do not interact 
with each other Jl|-^| . 



IV. PDF'S OF EIGENVALUE RATIOS IN TWO 
AND THREE DIMENSIONS 

We saw in the previous section that F(t, (), the PDF 
of the ratios Xi/g 1 ^, would remain the same in both Eu- 
lerian and Lagrangian frames. In this section we anal- 
yse the equations for these PDF's in two- and three- 
dimensional cases. Having in mind numerical simula- 
tions, we will write these equations using d — 1 indepen- 
dent variables. In the general case such reduction is done 
in Appendix 

Let us start with the two-dimensional case. It is now 
convenient to integrate the (^-function in (|l^) and work 
with the logarithm of the eigenvalue ratio as a new vari- 
able: x — ilog(Ai/A 2 ) = |(Ci — Ca)- The equation 
for F(t, x) then becomes 



d t F = (1 



F" + 

xx tanh(a;) 



F' 



(18) 



As expected, the rhs of Eq. (JiJ) becomes a full deriva- 
tive when multiplied by the Jacobian J{x) = 2sinh(ir). 
Note that the differential operator in the right-hand 
side of Eq. (|TJ) becomes a Legendre operator under 
the change of variables x = cosh(ir). This property 
is a consequence of integrability of the initial Hamilto- 
man H s (Eq. (§)), and will be of use m Sec. when we 
calculate the moments of passive vectors. 

The nature of the solution can be easily under- 
stood if we first consider only the advective term 
F' x l tanh(x). The characteristic of Eq. (18) then satis- 
fies x = l/tanh(x), which implies that F is advected 



to regions where 3> 1, and, for t — > oo, the asymp- 
totic solution can be found from ( |is| ) by approximating 
tanh(x) 1. The asymptotic is log- normal as expected. 

Note that the reflection symmetry x — > —x of Eq. (|T^ ) 
is just a consequence of the previously mentioned gen- 
eral symmetry Ai <-» A2, and does not add anything new. 
The function F must be initially chosen in such symmet- 
ric form. This is not so in the three-dimensional case 
that we now consider in more detail. 

In three dimensions, integrating the S- function 
in (HO) as before and introducing new variables, x = 
^\og[Xi/ A3) and y = ^log(A2/Aa), we obtain the equa- 
tion for F(t, x, y): 



F ' 



+ F" - 

; 1 xy 

sinh(x) 



F" 

yy 



+ 



d t F = (1 

.(-J— 

\tanh(a;) 

, 

\tanh(y) 2sinh(a;) sinh(y — x) 



2 sinh(y) sinh(x — y) 
sinh(y) 



F' 



F'„ 



(19) 



The normalization Jacobian for this PDF is J(x, y) 



sinh(x) sinh(y) sinh(x — y). 

The symmetry with respect to all permutations of 
eigenvalues Ai, A2, A3, leads to the following two symme- 
tries of the solutions of Eq. ([l9|): 

x —>■ —x, y — > y — x; and x «-» y. (20) 

Eq. (|l9]) posesses another (reflection) symmetry as well: 



x 



y 



-y- 



(21) 



which corresponds to the inversion of A1/A3, A2/A3, and 
does not follow from (|20|). Therefore a general initial 
distribution should contain both symmetric, F", and an- 
tisymmetric, F a , parts with respect to this reflection. 
The symmetries (^o| ) act as reflections (|2l]) on the points 
of the plane located on the lines y = 2x, y — x/2, 
and y = —x; hence the antisymmetric part of the PDF 
F a must vanish on these lines. 

Characterictic trajectories of Eq. (|l9|) are presented 
in Fig. |l|. The lines y = ±x, y — 2x, y — x/2, 
x = 0, and y = are combined in groups that are trans- 
formed by the symmetries (^0|) independently. Those 
groups correspond to sheet, tube, and strip volume de- 
formations as shown. Let us concentrate our attention 
on the sector x > 0, y < 0. Due to the symmetries ( |20| ) 
and (HJ), this allows us to understand the behavior of 
the PDF in the entire plane (x,y). Considering the 
characteristic trajectories (they advect F towards the 
line y = —x from both sides), or the flux of the conserved 
function F(x,y)\J(x,y)\ (calculated on the line y — —x, 
it is found to be directed from the semisector with pos- 
itive F a to that with negative F a ), one can show that 
the antisymmetric part of the PDF decays with time. 
The symmetry of the solution with respect to the sheet 
and tube configurations thus emerges asymptotically as 
t — * 00. 
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sheets 



FIG. 1. Characteristic trajectories of Eq. (|L9[). 
Solid lines correspond to sheet configurations, dashed 
lines to tubes, and dotted lines to strips. 

Numerical simulations performed for various initial dis- 
tributions concentrated in the region \x\ < 1, \y\ < 1, 
confirm that the PDF becomes symmetrized very fast, 
at times t ~ 1. In the region \x\ ^> 1, |y| ^ 1, far from 
the lines y — x, y — 0, and x — 0, the long-time (t 1) 
asymptotic is log-normal. This asymptotic can be easily 
obtained from Eq. ((L9|). 



V. PASSIVE VECTORS 



In this section we apply the developed formalism to 
passively advected vectors. Consider the evolution of the 
coordinates of a particle advected by the fluid: x % = 
x l (t,y), where y % is the initial position of the particle, 
i. e. x*(0,y) = y l . An infinitesimal contravariant vec- 
tor a 1 changes under such coordinate transformation as 
follows: a l (t, x) = (dx t /dy k ) a§(y)- In order to find the 
mean of any object constructed out of a 1 , we have to av- 
erage it with respect to the initial distribution of a (y) 
and with respect to all realizations of the random velocity 
field £ k . The latter averaging can be done via the PDF's 
for co- and contravariant (metric) tensors. Let us assume 
that the initial distribution of the vector a 1 is Gaussian, 
isotropic, and independent of y: (a^a^ = . As an 
example, consider the moments A n = (|a| 2 ): 



A n = ((a • g ■ a ) n >, 



(22) 



where g is the contravariant tensor advected by the fluid 
and with the initial condition g^(0,y) — The distri- 
bution of this tensor can be found in the same way as that 
of the covariant tensor, and is discussed in Appendix |^. 
Moments of a covariant vector ai can be found using ex- 
actly the same formula (p2|), with g now the covariant 
tensor. 

To simplify the formula (|2^), we note that the eigenval- 
ues of the matrix g can be expressed as A; = g x / d exp(Ci). 
Therefore, for all n, g~ n l d Tr (g n ) depend only on Q, and 
are independent of the determinant g. Since the initial 
distribution oft; is S(g) = 5(g—l), we can average powers 
of g independently and obtain: 



A n = f d (n,t)(g n ? d ), 
f d (n,t) = (g- n ' d (ao-S-ao)"}, 



(23) 



where the functions fd(n, t) do not depend on the statis- 
tics of the determinant, and are, therefore, universal. 
These functions are the same in the co- and contravariant 
cases, and in both Eulerian and Lagrangian frames. The 
only parts of the moments A n that are non-universal are 
the averages of the determinant. These averages can be 
calculated exactly using formulas (|l|), @, and (|c|): 



(9%° = (.9 



(,9 s 



(fJ 



s^contra _ s (2s-l)ft 
L ' 
s\contra _ s(2s+l)ft 
E — L 



(24) 
(25) 



The universal functions fd{n,t) can, in fact, be easily 
calculated directly (cf . P,p yT^ , p^ | ) , if one starts from the 
equation for advection of a passive vector a l (t,x): 



d t a i + ea: k -e k a k 



0, 



(26) 



where statistics of £ l (t, x) are given by (g). However, for 
methodical purposes, we prefer to rederive this result us- 
ing the technique of Jack polynomials. While it is also 
quite simple, it illustrates the general method that can 
be applied to finding moments of any passively advected 
tensor. At the end of this section, we show, e. g., how 
moments of a bilinear form a l b k can be caclulated. 

Formula (^3|) can be further simplified if we do the av- 
erage with respect to the distribution of <Zg. Introducing 
the generating function 

Z{p) = (exp {(3 g- x / d (a • g ■ a ) } ) , (27) 

we represent fd as follows: 

~d n Z((3) 



fd(n,t) 



d(3 n 



(28) 



0=0 



The Gaussian average with respect to the initial distri- 
bution of the vector can now be easily done, resulting in 



-1/2' 



(29) 



where the remaining averaging is with respect to the 
statistics of Ci- The PDF of the C's is F(Q\ J(C)|<HE CO 
with the initial condition <5(Ci) • • • S((^d)- It is important 
that the function that is being averaged in ( p9| ) is the 
generating function for a particular class of Jack polyno- 
mials, that are eigenfunctions of the self-adjoint Suther- 
land operator Hs in (||). Therefore, all functions ( |2^ ) 
can be found exactly in the general case. The appropri- 
ate calculation is carried out in Appendix The answer 
is: 



fd(n,t) 



exp 



f 



n(2n + d)(l + a)t 



(30) 



5 



where we denote: (d/2) n = (d/2)(d/2+l) • • • (d/2+n-l). 

In the two-dimensional case the corresponding result 
can be obtained in a rather simple manner, which never- 
theless illustrates the main idea of the general derivation. 
In order to do this, we notice that the generating func- 
tion Z(f3), expressed in the two-dimensional case in terms 



of x = i(Ci — C2) (see Sec. IV), coincides with the gener- 



ating function for the Legendre polynomials P„(cosh(a;)), 
and, therefore, f2(n,t) = n! (P„(cosh(x))). The average 
can now be completed with the aid of Eq. (|l8|). Mul- 
tiplying it by \J(x)\ P n (cosh(x)), integrating by parts 
twice, and using the equation for the Legendre polyno- 
mials, P"(p)(p 2 - 1) + 2^0) = n(n + l)P„(/x), we get: 

f2(n, t) = (1 + a)n(n + 1) /^(n, t), whence 
f 2 (n,t) = n!exp{n(n+ 1)(1 + a)t}, (31) 



which is in agreement with (pO[). 

As an example, consider moments of a magnetic field 
advected by the fluid. The contravariant vector in this 
case is B l /p, where p is the density of the fluid. Let us de- 
note the moments of B l as H n = (|B| 2 ™). Recalling that, 
in the contravariant case, g — l/p 2 , we get from (|23|): 



fd(n,t)(g 



-n{d-l)/dy ontl '& 



(32) 



where for the g average we use the formula (|25|) in Eule- 
rian frame, or (|24|) in Lagrangian frame. 

An analogous derivation can be carried out for a co- 
variant vector, e.g., gradient of a passive scalar V9. For 
its moments C n = (|V#| 2 ™), we find: 



C n = h{n,t)(g n l d Y 



(33) 



where for the g average we use formulas (24) or (|25|) de- 
pending on the frame of reference. 

On passively advected tensors 

We now briefly demonstrate how one can calculate ex- 
actly the moments of a passively advected higher-rank 
tensor T. Suppose that we are interested in some mo- 
ment (T m ). After averaging with respect to the ini- 
tial distribution of T, we are left with a combination 
of Tr fc (<7"), which are polynomials of degree nk in the 
eigenvalues of the metric g. But any symmetric polyno- 
mial of degree m can be expanded in Jack polynomials 
of degree m, which can then be averaged exactly. The 
result will therefore be a linear com bination of exponents 
growing at the rates given by flB10| ) . 

For example, consider a contravariant bilinear 
form a l b k , where a 1 and b k are initially independent 
Gaussian random vectors, (aga§) — (&o^o) = <^ fe > an d 
find its second moment B2 = ((a-b) 2 ) = ^Tr((p)). 
Then 



Bo 



1 d 



2 ) = (9 Vd ) 



(J(.2,0)) - o ( J (l,l)) 



(34) 



where polynomials J(2,a) an d J(i,i) are constructed 



' (2) 

in (B4). The corresponding eigenvalues are E y ' 



(d + 4)(d - l)/d and 
from (BIO). The answer is: 



2,0) 

(d 2 - A)/d, as follows 




(35) 
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APPENDIX A: PDF OF EIGENVALUE RATIOS 

The ^-function in ( |Io| ) can be integrated over, 
and Qi,...,Qd reduced to d — 1 independent variables, 
viz. the logarithms of the eigenvalue ratios: x n = 
\ log(A n /A<z) = 5(Cn ~ (d)- In these variables, the equa- 
tion for F becomes: 



d t F=(l + a)\J2 



d 2 F 



E 



1 



OF 



z - ' dx 2 n ^ tanh(a;„) dx n 



d-1 

4e 



d 2 F 



2 Q^CfidXfi 



1 d-1 



4 sinh(x„ - x rn ) 



sinh(x m ) dF 



sinh(a; n ) dx r , 



sinh(x„) dF 
sinh(a; TO ) dx n 

(Al) 



The last two terms correspond to interactions between 
different a;'s and only enter for d > 3. The normalization 
rule now is: 



2(d+2)(d-l)/2 



d-i 



n i sin M 



x Y[ I smh(x n )\dx n F = 1. 



(A2) 



71=1 



This form of the equation for F is most convenient 
for numerical solution and for geometric analysis such as 
that of Sec. Ev. 



G 



APPENDIX B: JACK POLYNOMIALS 

Jack polynomials J^{x\, . . . ,Xd', ol) of degree m are ho- 
mogeneous (of degree m) polynomials, depending on d 
variables x%, and symmetric under all permutations of x^. 
They depend on a parameter a and are labeled by par- 
titions [i of an integer number m. 

Partition fi of m is a non-increasing sequence of in- 
tegers: /i = (fii > ... > fid) € Z> , such that m = 
fix + . . . + fj,^. The polynomials J^(a;;a) vanish if the 
number of parts ?(/z) is greater than the number of vari- 
ables d. Consider two partitions fi and A of the same 
length = 1(A) = d. One writes that A > /i 

if Ai + . . . + A, > fix + . . . + fii for each i < d. This 
defines the so-called natural (or dominance) ordering on 
partitions. 

In order to give a formal definition of the Jack polino- 
mials, first define the monomial symmetric function m^, 
corresponding to the partition fi: 



= E< 



(Bl) 



where the summation is over all permutations 
of fix, . . . , Hd- The Jack polynomials must, by defini- 
tion, be represented as: 



J\(x;a) = ^2 ux^m^, 



(B2) 



and be the eigenfunctions of the Calogero-Suthcrland 
Hamiltonian: 



H 



(a) 

2 — 1 



d \ 2 2 * x? d 



Ehf- + 1 E — f- 

^ — ' \ rlfr.j I rv — < Xi — I, OXi 



(B3) 



All coefficients can be found recursively in terms 
of u\\ with the aid of this definition [ji"7f . 

Let us use this definition to construct the Jack poly- 
nomials for m — 2 and a — 2. The corresponding par- 
titions are (2. 0, 0, . . .) and (1, 1,0,.. .). Using the first 
condition (B2), we write: 



J(2,o) (x; 2) = ^2 x i + A E XiX i ' 



i=i 
a 



J(ia)(x;2) =E a; 



(B4) 



i<j 



The coefficient A must be found from the requ irement 
that the polynomials be eigenfunctions of (B3), which 
gives A = 2/3. 

The eigenvalues (energies) corresponding to Jack poly- 
nomials are: 



d 2 d 



(B5) 



The energies (B5) depend on particular partitions fi. Any 
symmetric polynomial of degree m can be expanded in 
Jack polynomials of the same degree m. Of all the other 
properties of the Jack polynomials we will need the fol- 
lowing: 



i 



(l-xtyj) 1 /* 



22b fJ ,(a)Ji l (x;a)J tl (y;a), (B6) 



where the summation is performed over all possible par- 
titions fi of all non- negative integers, and b u (a) are some 
expansion coefficients that can be foun d in J17| . For our 
purposes we will need the formula (B6) with the set {yj} 
consisting of only one variable. In this case the expansion 
takes the form: 



d ^ x 

n ( i.^i/ a = E^w(';°). 



(B7) 



m=0 



where fi = (m) is a partition consisting of only one ele- 
ment. Q( m )(x; a) stand for the properly normalized Jack 
polynomials. The explicit expression for Q^ m ^(x;a) is as 
follows: 



Q(m)( x l a ) = E 



Kii...<j 



Xi-, Xi. 



■x im , (B8) 



where 9 = 1/a, q\ = JJ={n\i n = 1} is the multiplicity 
with which the number I = 1,2, ... ,d appears in ix . . . i m , 
and [9) q = 6(9 + l)---(9 + q-l). 

To use these results we need to transform our Hamil- 
tonian (in the square brackets in (^)) to the form (|B3|). 
Changing variables to A, = exp(£,) = Xig^ 1 / 4 , we get: 



He = - 



1 d 

zE 



d^dQdQ 



+ 



1 d 

It, 



tanhi(Ci-Ci) dQ 



d 2 

d 



d 



H 



(2) 




-E( s 



_d_ 



(B9) 



For any Jack polynomials of degree m, the corresponding 



eigenvalues E^ of the Hamiltonian (BE) are: 



m 
~d 



(d-1) 



■ m. 



In particular, the energy of Q( m )(A; 2) is 
d- 



1 / d 
— mim+- 



(B10) 



(Bll) 



;.=i 



We now notice that the averaged Jack polynomi- 
als <3(„)(A; 2) and the functions fd(n,t)/n\ have the same 
generating function (see (p9h and (B7)), whence 



7 



f d (n,t)=n!(Q {n) (\;2)). 



(B12) 



Since (Eq. ®) d t F = 2(1 + a)HsF, where # s is self- 
adjoint with respect to the measure (|l^), fd(n,t) satis- 
fies: fd(n,t) — 2(1 + a)E n fd(n,t), the solution of which 
(with correct initial condition) is the expression (|3C 



APPENDIX C: PDF FOR CONTRAVARIANT 
TENSOR 

The derivation of the main equations for the case of 
a contravariant tensor is quite similar to the case of the 
covariant tensor. Here we just explain the origin of the 
difference and write out the main results. The dynamical 
equation in the contravariant case reads: 

dtv> ij + - e k v hj - e k ^ ik = o. 



(Cl) 



The derivation of the master equation can be carried out 
the same way as in the covariant case and results in dif- 
ferent coefficients in the g part of Eq. (^|). The £ part 
remains intact. This is not surprising, since the tran- 
sition from tpij to ip lJ does not change the ratios of 
the eigenvalues, but results only in the inversion of the 
determinant: g — > 1/g. Accordingly, the equation for 
the PDF of a contravariant tensor, P(g,(), can be ob- 
tained from Eq. (||) by substituting P = Pg d+1 , g = 1/g. 
The resulting equation for the PDF of the determi- 
nant, S(g), is (dropping the overtildes): 



d t S = 2g 2 ^ + (2d 



3) 5 — + -d(d+l)S, (C2) 



where time has been rescaled by the factor of 7 as 
in Sec. 



Ill 



When using this equation, we should re- 
member that g now satisfies the same equation a s 1/ p 2 , 
where p is the density of the medium. Eq. (C2) is 
written in the Eulerian frame. For completeness, we 



write down the solution of Eq. (C2) with initial distri- 
bution S{0,g) = S(g-1); 



S(t,g) 
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-(d+l)/2 



■ exp 



(log(g) -ity 

8jt 



(C3) 



The Lagrangian analogue of (C3) is obtained via multi- 
plication by p — 1/ ^fg. 
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